Stringent constraints on cosmological neutrino-antineutrino asymmetries 
from synchronized flavor transformation 
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We assess a mechanism which can transform neutrino-antineutrino asymmetries between flavors 
in the early universe, and confirm that such transformation is unavoidable in the near bi-maximal 
framework emerging for the neutrino mixing matrix. We show that the process is a standard 
Mikheyev-Smirnov-Wolfenstein flavor transformation dictated by a synchronization of momentum 
states. We also show that flavor "equilibration" is a special feature of maximal mixing, and carefully 
examine new constraints placed on neutrino asymmetries. In particular, the big bang nucleosyn- 
thesis limit on electron neutrino degeneracy £ e | < 0.04 does not apply directly to all flavors, yet 
confirmation of the large- mixing- angle solution to the solar neutrino problem will eliminate the 
possibility of degenerate big bang nucleosynthesis. 
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I. INTRODUCTION 

As is well known, the observational successes of big 
bang nucleosynthesis (BBN) are one of the pillars of stan- 
dard cosmology 0. If one assumes three standard neu- 
trino flavors, then the only free parameter is the baryon 
to photon ratio ng/)i 7 . The value ns/n 7 ~ 5 x 10~ 10 
predicts light-element yields of 2 H, 4 He, and 7 Li that 
are in excellent agreement with observations. As is often 
noted, this is particularly remarkable because the abso- 
lute yields of these elements span several orders of mag- 
nitude. This consistency implies that the post-BBN pro- 
cessing of the light elements is largely understood, and 
that one does not require new aspects of particle physics 
beyond the standard model (SM) that would materially 
affect BBN. 

The basic consistency of our picture of the early uni- 
verse is even more impressive when one considers other 
recent cosmological measurements. Observations of the 
acoustic peaks in the angular power spectrum of the cos- 
mic microwave background (CMB) give strong evidence 
that Vltotai = 1-04 ± 0.06 g, i.e., the universe is flat, 
as predicted by inflation. Taken together with measure- 
ments of clusters of galaxies and the high-redshift type-la 
supernovae (SNIa) data, a mutually consistent picture || 
with flmatter — 0.3 and fliambda — 0.7 is obtained. Addi- 
tionally, the CMB data indicate that Qbaryon — 0.04, in 
excellent agreement with the BBN observations. Recent 
data on the clustering of galaxies also yield consistent 
values of ^matter and flbaryon @ • This agreement of the 
combined data is all the more impressive because bary- 
onic matter is such a small fraction of the total energy 
density of the universe, and because the BBN and CMB 
data reflect measurements of Qbaryon at very different 
epochs (~ 10 2 s and ~ 10 13 s after the big bang, respec- 
tively) . 

Nevertheless, from a particle-physics point of view, 
these impressively measured quantities are all totally un- 
explained, both in their values and their nature (e.g., 



though we know that the particle dark matter is not 
part of the SM, we do not know what it is). Just as 
in accelerator-based particle physics, the underlying be- 
lief is that more and more precise measurements will lead 
us to the necessary clues on how to generalize the SM. In 
particular, the baryon-antibaryon asymmetry of 5 x 10~ 10 
remains a mystery, and is certainly an important clue for 
understanding the universe at temperatures at least as 
high as the electroweak scale. 

Naturally, attention is also focused on the lepton- 
antilcpton asymmetry of the universe. General consid- 
erations indicate that B — L may be conserved, so that 
the lepton asymmetry is nr, /n 7 ~ 5 x 10 -10 as well. How- 
ever, there are certainly viable models in which the lep- 
ton asymmetry can be much larger and if confirmed, 
would be a very important clue. Given constraints on 
charge asymmetry, any large lepton asymmetry would 
have to be hidden in the neutrino sector. Though the 
baryon asymmetry can generically be limited to be less 
than 10~ 8 simply to not overdose the universe, no similar 
constraints exist in the lepton sector for light neutrinos. 

Since neutrinos and antincutrinos should be in chem- 
ical equilibrium until they decouple at a temperature 
T ~ 2 MeV, they may be well described by Fermi-Dirac 
distributions with equal and opposite chemical poten- 
tials: 
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1 + exp(p/T - £) 



(1.1) 



where p denotes the neutrino momentum, T the tem- 
perature and £ is the chemical potential in units of T. 
(There is a tiny non-thermal perturbation that occurs at 
the epoch of e + e~ annihilation at T ~ 0.3 MeV, which 
we can ignore.) The lepton asymmetry L a for a given 
flavor v a is related to the chemical potential by 



12C(3) 



so 



(1.2) 
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where £(3) ~ 1.202. Even enormous values of £ Q ~ 1 have 
been allowed observationally. This is so distant from the 
naive SM prediction that any measured nonzero value 
would be very important. Interest in searches for such 
large values of £ Q is also driven by the fact that we ev- 
idently have much to learn about the neutrino sector. 
Previously, most attention was devoted to the possibil- 
ity of significant mixing of the SM active neutrinos with 
light sterile neutrinos, which can generate large lepton 
numbers L ~ 1 |(|. 

A general approach to setting limits on £ Q arises be- 
cause for very large degeneracy, the effective number of 
neutrinos is increased from the standard model predic- 
tion by 



(1.3) 



This increases the expansion rate of the universe, chang- 
ing the CMB results by magnifying the amplitude of the 
acoustic peaks. For all flavors, the bound |£ a | < 3 has 
been obtained from the CMB alone Jjj. Note that the 
sign of £ a is unconstrained. With future CMB data, 
these limits may be reduced to |£ a | < 0.25 or less A 
much stronger limit can be placed on £ e with |£ e | < 0.04 
because of its effect on setting the neutron to proton ra- 
tio prior to BBN by altering beta equilibrium. 1 If at 
the same time £ MiT are large, this effect can be partially 
undone by the increased expansion rate, leading to the 
often-quoted bounds [|[ [L0| 



0.01 < i e < 0.22, 

l^.rl < 2.6, 



(1.4) 
(1.5) 



where the upper limits are obtained only in tandem. 

There are now three types of evidence for neutrino os- 
cillations: solar neutrino [ pT| v e — > v^, v T with large (but 
not maximal) mixing angle and 8m 2 ~ 10~ 5 eV 2 , at- 
mospheric |l2| ] neutrino , — * v T ,v T with maximal 
mixing and 8m 2 ~ 10~ 3 eV 2 , and the Liquid Scintillator 
Neutrino Detector (LSND) jl3| neutrino — » v e with 
a very small angle and 8m 2 ~ 1 eV 2 . It is not possible 
to accomodate these three signals with only three neu- 
trinos, as there are only two independent mass-squared 
differences. A possible fourth (sterile) neutrino can be 
invoked to create a new 8m 2 , but now that the solar 
and atmospheric neutrino data indicate the appearance 
of active neutrino flavors, there is a problem of where 
to incorporate the required mixing with the sterile neu- 
trino. While four-neutrino models may still work, it is 
only with difficulty, both in fitting the oscillation data 
(see, e.g., Rcf. |Q) and through the effects on BBN (see, 
e.g., Ref. p5|). The LSND signal will be conclusively 



confirmed or refuted by the MiniBooNE experiment [[16| . 
For simplicity, we consider just three active neutrinos, 
and neglect the LSND result (of course, if it is confirmed, 
a major revision will be necessary). 

In such a three-neutrino framework, Lunardini and 
Smirnov |rj| suggested that the large mixing angles im- 
plied by the present data may transfer any large asymme- 
try hidden in ^ T to £ e well before the beta-equilibrium 
freezeout at T ~ 1 MeV (see also Savage, Malaney and 
Fuller 18 ). Thus, the stringent BBN limit on £ e might 



1 Beta equilibrium is between the weak interactions n+v e <-» p+e~ 
and p + v e <-» n + e + . Positive §e increases the v e abundance 
relative to u e , forcing equilibrium towards lower n/p. 



apply to all three flavors, improving the bounds on £ MjT 
by nearly two orders of magnitude. 

This proposal was recently studied in detail by Dol- 
gov, Hansen, Pastor, Petcov, Raffelt, and Semikoz (DH- 
PPRS) pS| ]. They found that close to complete trans- 
formation of asymmetries £ M and £ r to £ e was obtained. 
This is an important result, as it excludes the possibil- 
ity of degenerate BBN j2C|, and is the strongest limit on 
the total lepton number of the universe and is likely to 
remain so for the foreseeable future. In this article we ex- 
amine the DHPPRS result, show it as the result of a syn- 
chronized Mikheyev-Smirnov-Wolfenstein (MSW) trans- 
formation and establish its robustness through physical 
and numerical insight into the dynamics. We assess how 
the results depend on the input parameters and consider 
more exotic physical scenarios that might affect the re- 
sults. 



II. TWO-FLAVOR DENSITY MATRIX 
EQUATIONS AND SYNCHRONIZED MSW 

In this section we consider a mixed neutrino statisti- 
cal ensemble in the early universe, with initial ncutrino- 
antincutrino asymmetries which are not equal among fla- 
vors. We show that this ensemble behaves as a syn- 
chronized system following a single effective momentum 
state that undergoes an MSW transformation given large 
mixing angles. In describing neutrino flavor evolution 
in dense environments such as the early universe, one 
must use a density matrix description if the neutrino self- 
potential is large or if decohering collisional processes are 
significant. Where collisional processes are not impor- 
tant, as is the case for some examples we shall consider 
here, the evolution is coherent. A useful parametrization 
of the density matrix equations is the Bloch form [2l| . 

In an environment such as the early universe, where the 
potential arising from neutrino-neutrino forward scat- 
tering (the neutrino self-potential) is important, active- 
active mixing is substantially different from active- 
sterile mixing. Forward scattering processes of the type 
v aip) ~ * v p(p) l ea d to refractive index terms which are 
off-diagonal in the flavor basis {a, (3} |p2[ . A useful and 
interesting casting of the Bloch formalism for pure active 
neutrino mixing was done by Pastor, Raffelt and Semikoz 
p3j , which allows an interpretation via analogy with the 
precession of coupled magnetic dipoles. The analysis of 
Ref. j23| considered the case of constant density, in the 
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absence of a background medium other than that pro- 
vided by the neutrinos themselves. We shall have need to 
extend this description to include a background medium 
of charged leptons of a density that varies with time (or 
temperature). One particularly interesting feature first 
revealed clearly in Ref. pjjfl , is that the neutrino self- 
potential (that is, the potential due to neutrino-neutrino 
forward scattering) does not, in general, suppress flavor 
oscillations, as one might have naively expected by anal- 
ogy with the potential from, say, a background of charged 
leptons. This is in stark contrast to the case of active- 
sterile oscillations, where the effect of an asymmetry be- 
tween the active neutrinos and antineutrinos is always 
to suppress mixing angles. Even for relatively small de- 
generacies, the asymmetry term dominates the evolution 
and thus delays transformation of such asymmetry from 
the active to the sterile flavor. 

For active-active oscillations, no such simple mixing 
angle suppression occurs as the neutrino asymmetry en- 
ters both the diagonal and off-diagonal terms in the ef- 
fective Hamiltonian. These terms have the effect of syn- 
chronizing the ensemble, 2 resulting in collective behavior 
resembling the evolution of a single momentum state in 
the absence of the self-potential. The mixing angle for 
this collective oscillation is determined essentially by the 
background medium of thermal charged leptons. When 
the density of this background decreases with tempera- 
ture, the neutrinos evolve adiabatically from their initial 
flavor into vacuum mass eigenstates. For large-angle mix- 
ing, this implies significant flavor transformation. 



A. Formulation 

We express the mixing between two neutrino mass and 
flavor eigenstates as 



JdV(27r) 3 /fe0) = 3C(3)T 3 /4^ 



cos Vq vi + sin Wo V2 i 
sin 6qi>i + cos #0^2, 



(2.1) 



where, in general, v* denotes some linear superposition 
of and v T , as we shall explain later. We parametrize 
the two-flavor neutrino density matrix in the form 

p(p) = ( P pH)=\ i P oiP) + ° ■ P (P)] - ( 2 - 2 ) 

and similarly for the antineutrinos, where we refer to 
P(p) as the neutrino "polarization" vector. These quan- 
tities are most usefully normalized such that 

P(p) mMal = [/ e (p,6) ~ U(P,^)} , (2-3) 



2 Note that this forward-scattering induced synchronization is un- 
related to the synchronization effect of Ref. ]24l which arises from 
rapid flavor-blind collisions. It is remarkable that both collisions 
and forward scattering lead to synchronization effects in the case 
of active-active oscillations. 



where 7ig q 

In the absence of collision terms, Pq(p) and the mag- 
nitude of P(p) remain constant and the full evolution 
equations for two mixed active flavors in the early uni- 
verse are 

d t P P = +A P x P p + a(3 - J) x P p , (2.4) 



c\P E 



a(J - J) x P 



V 1 



where P p denotes the polarization vectors for the neutri- 
nos of momentum p while J denotes the corresponding 
quantity integrated over momentum such that 



1 (2tt)3 ^ 



(2.5) 



Vectors with an overbar refer to the antineutrino quanti- 
ties throughout. With the normalization we use here, the 
length of the individual P p vectors and that of J do not 
redshift with temperature. The coefficient of the second 
term is a = v / 2Gi?"< cq . Time t and temperature T may be 
interconverted via the expression t ~ 1.15 s (T/MeV)~ 2 . 
Decohering collisions (damping) of the system at high 
temperatures forces the neutrinos into unmixed flavor 
states. We shall assume zero initial and a finite initial 
taken to have a negative sign. 3 Therefore, the initial 
alignment of the P p are along the +z axis, and P p are 
along the —z axis. 

Equations (2.4) are equivalent to the precession 
of magnetic dipoles in two "magnetic fields" : the 
momentum-dependent A p and the integrated neutrino 
self-potential a(J — J). The effects of A p , as we shall 
show, are straightforward, but the neutrino self-potential 
makes the system non-linear by explicitly coupling each 
momentum mode to the evolution of every other mo- 
mentum mode. An intuitive description of the evolution 



in Eq. (2.4) for a constant-density system without mat- 
ter effects was provided in Ref. |2l| . The issue of the 
synchronization of the system has also been studied in 
Ref. |!J. 

In general, the "magnetic field" vector A p includes 
contributions from vacuum mixing, a thermal potential 
from the charged-lepton background, and a potential due 
to asymmetries between the charged leptons 



A P = A P + [V T (p) + V 1 



(2.6) 



Vacuum mixing is incorporated by 



A p = ^^(sin2y x-cos26>oz), (2.7) 
Zp 

where Sitiq = to 2 . — m\ and 0$ are the vacuum oscillation 
parameters. 



3 For the opposite sign, one simply reverses the directions of the 
polarization vectors. 
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The thermal potential from finite-temperature modifi- 
cation of the neutrino mass due to the presence of ther- 
mally populated charged leptons in the plasma is 



V T {p) = 



8V2G F p 
3m^ 



+(E l+ )n l+ ) . (2.8) 



The neutrinos also contribute a thermal self-potential 
similar to the form of th e se lf-potential on the right- 
hand side (RHS) of Eqs. (0 @, but unlike a ~ G F , 
the thermal self-potential goes as Gp. Unless the initial 
asymmetry is of a size much too small to be interesting 
here, the G 2 F term is negligible by comparison with the 
order Gf self-potential and thus unimportant in deter- 
mining the dynamics of the system. 

The background potential arising due to asymme- 
tries in charged leptons is nonzero only for electron 
neutrinos to maintain charge neutrality of the baryon- 
contaminated plasma: 



B= {±V2G F (r 
10 



n e +) fort/ e ^^ IjT , 

for V^^Vr, 



(2.9) 



where +(— ) is for neutrinos (antineutrinos) . Due to the 
smallness of the baryon asymmetry relative to number 
densities of thermalized species, this term is always neg- 
ligibly small relative to the vacuum vector A p and the 
thermal potential V T , and so we may take A p = A p . 

In the absence of the neutrino self-potential it is pos- 
sible to define 



A — 
zip — 



2p 



-(sin26> m x — cos26> m z), 



(2.10) 



where <5m^ and m are the matter-affected oscillation 
parameters. When the self-term must be included, the 
nonlinearity of the problem makes the notion of a matter- 
affected mixing angle more subtle. 



B. Synchronization 



Taking the difference of the Eqs. (2T), integrating over 
momenta, and defining a collective polarization vector 
I = J — J one obtains 

d t I ~ A cff x I, (2.11) 

where the appropriate effective "magnetic field" is 



A cff ~^ / A,, (P., + P,,i-I. 



(2.12) 



In fact, Eqs. ( 2.11 ), ( 2.12| ) are exact only when I | 
(Pp + Pp). The vector I thus precesses slowly about A c ff . 
Since the self-potential dominates Eqs. (2.4), the individ- 
ual polarization vectors P p all precess rapidly about I, 
and, if initially aligned, are held together. The various 
vectors are illustrated in Fig. |l|. 
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FIG. 1: Vector precession diagram. The angles and magni- 
tudes of the vectors are not to scale but have been exaggerated 
for clarity. For the situation of interest, the magnitude of I is 
much greater than that of A c ff . When this condition holds, it 
is a good approximation to describe the evolution of the po- 
larization vectors for the individual momentum modes P p as 
a precessing about I. The vector I then precesses about A e ff, 
in the manner of a single momentum mode in the absence of 
the self-term. For asymmetries between neutrino flavors in 
the early universe, A e ff and I are both initially aligned with 
the z axis, and, for maximal mixing, both adiabatically evolve 
to align with the x axis. 



We assume for simplicity that the initial asymmetry 
resides in a single flavor. Although the coefficient of the 
self-term makes it dominant, the flavor transformation is 
actually determined by the evolution of A c ff. In fact, if 
one leaves out the self-term altogether, the synchroniza- 
tion is of course lost, but the average flavor evolution of 
the system is almost completely unchanged. 

Let us first consider the sim ple linearized case where 
the self-potentials in Eqs. (2.4) vanish. In this case, the 
polarization vector of each momentum mode P p precesses 
about its own respective A p : 



d t P P 



(2.13) 



c^Pp — Ap x Pp . 



If one follows only the average momentum (p/T) ~ 
3.15, Eqs. ( |2. 13 ) are simply two linear equations with 



a straightforward solution. Recall that each P p is ini- 
tially aligned along the +z axis. At high temperatures, 



\V L | > |A p |, and thus from Eq. fl2.6| ) the A p point in 
the — z direction. As the temperature of the universe 
decreases, \V T \ ~ T 5 decreases and the vectors A p will 
slowly rotate from the — z direction toward the +x. di- 
rection and the angle that A p subtends with the z axis 
will asymptote to 26*o- This effect is a straightforward 
MSW transformation of the asymmetry: the initial neu- 
trino number excess in one flavor evolves from a mass 
eigenstate in matter (modified by the thermal potential) 
to a vacuum mass eigenstate with different flavor con- 
tent. And, since Eqs. ( p,13j ) are decoupled, the average- 
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momentum mode will describe the collective evolution of 
the entire system. 

In the substantially more involved system, including 
the self-potential ( |2.4| ), each momentum mode is coupled 
to all other momenta through the self-term. Therefore, a 
simplifying average-momentum technique is poorly jus- 
tified. However, if one blindly drives f orw ard with an 
average-momentum evolution of Eqs. (2.4), one luck- 
ily recovers (nearly) the correct behavior. For the full- 
momentum case, including the self-potential, the I vector 
will also initially be al igned with the z axis. Using the ap- 
proximate Eqs. (2.11) and making the as-yet unjustified 
assumption that A e g follows the average of the A p , the 
synchronized system will undergo the MSW transforma- 
tion at the exact same temperature as the over simplified 
case (2.13). This is what we found numerically. 



The fact that the neutrino flavor system evolves to the 
same end-state with and without the self-term appears at 
first absurd. However, one must keep in mind that the 
neutrino self-potential in a purely active neutrino sys- 
tem affects the evolution drastically differently than the 
familiar flavor-diagonal potentials in the matter Hamil- 
tonian present, for example, in evaluating the solar MSW 
effect and active-sterile mixing in the early universe. The 
neutrino self-potential, as the dominant precession term 
for each momentum, only plays the role of forcing each 
momentum mode to follow the collective vector I. 



To explore the behavior of the system and verify the 
appro xi matio ns in arriving at the collectiv e e quations 
(2.11), (2.12), we numerically integrate Eqs. (2.4), which, 
again, explicitly couple the full thermal distribution of 
momenta to the quantum mechanical evolution of each 
momentum mode. Because of the drastically different 
time scales over which the terms A p , V T (p) and a(J — J) 
evolve, the system is a set of stiff nonlinear differential 
equations and therefore requires careful treatment. 



In Fig. ||a we show the evolution with the full equa- 
tions ( |2.4| ) of a representative two-flavor v e and system 
with the best-fit solar large mixing angle (LMA) param- 
eters. Shown is the angle between each A p and z, which 
would determine evolution of each P p in the absence of 
the self-term. The actual angle between P p and z is 
shown in Fig. ^|b, displaying the stunning synchroniza- 
tion of all momentum modes, to roughly the orientation 
of A p at the average momentum (p/T) ~ 3. Figure || 
shows the tiny magnitude of the angle between P p and 
I, which is the result of the synchronization. 

Now, we justify why taking the evolution of A c ff to 
be effectively that of the average of A p luckily provides 
the nearly the correct evolution. For the case where col- 
lisional damping may be neglected and assuming that 
both synchronization and the close alignment of P p + P p 
with I holds, we may explicitly calculate A c ff which de- 
scribes the MSW-like transition (the expressions for A c g- 
and Psync were independently calculated in Ref. (p7f). We 
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FIG. 2: The angle between A p and the z axis is shown in 
the upper panel, in color, as a function of the temperature of 
the universe (horizontally) and across the neutrino spectrum 
(vertically). In the lower panel, the angle between P p and 
the z axis is displayed in the same fashion. As detailed in the 
text, all P p ignore the momentum dependence of A p and are 
dramatically synchronized to a single effective momentum, 
Psync/T ~ 7r. That is, all of the P p follow the orientation 
(i.e., have the same color) of A p at p/T ~ it, shown with a 
white horizontal dashed line. 



find 

A cf f ~ f- 



3/2 



2T \n 2 +^ 



sin 26» x + (- cos 26» + Z) z , 
(2.14) 



where 



Z 



2T 



SrriQ \p/T J 



(2.15) 



Note that Z is negative so th ere is no resonance. It is 
helpful to reexpress Eq. ( 2.14 ) in terms of effective "syn- 
chronized" oscillation parameters as 



A e g — A S y nc (sin 2$ S y nc x cos 2$ S y nc 2 



(2.16) 
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FIG. 3: We show the angle between individual P p and I as 
a function of temperature of the universe (horizontally) and 
the neutrino spectrum (vertically). The angles are extremely 
small, indicating the degree of synchronization. 



where the synchronized oscillation frequency is given by 
jm% ( 3/2 



e- 



2T \n 2 + e 



sin 2 20 o + (-cos26» + Z) 2 . 

(2.17) 



The size of A sync , which is proportional to an overall 
factor of £, is not important in determining when the 
MSW-like transformation takes place. 4 It is the mixing 
angle 8 sync that is important in describing the transfor- 
mation resulting from the evolution into vacuum mass 
eigenstates. This angle (i.e., half of the angle between 
A e ff with the z axis) is given by 



sin 2 26. 



sin 2 26>n 



sync 



sin^ 26> + (-cos2(9o + Z) 2 



(2.18) 



and thus we find it is the mixing angle which would cor- 
respond to the momentum state 



Psync 



HyJ 1 + £ 2 /27T 2 ~ 7T. 



(2.19) 



This is one of our principal results, and indicates a re- 
markable coincidence. Namely, the apparently identical 
evolution for the synchronized system including the self- 
term and that found with a vanishing self-term only re- 



sults from the fact that average momentum for a rela- 
tivistic Fermi gas (with small chemical potential) 



(P/T) 



7tt 4 



180<(3) 



3.15 



(2.20) 



is approximately 7 r, the effective momentum of the syn- 
chronized system (2.19) with £ -C tt. One can observe in 
Fig. H that the effective mixing angle (the angle of A p 
with respect to the z axis) does in fact correspond to 
the way in which the state p/T ~ 3 would evolve in the 
absence of the self-term. 

It is clear that 9 sync depends only very weakly on the 
size of the initial asymmetry, and in particular the trans- 
formation will occur at almost the same temperature for 
any plausible initial £. Additionally, if £ were very large, 
even a very small degree of flavor transformation would 
be sufficient to upset successful BBN. We have also nu- 
merically integrated the system for several initial asym- 
metries, and verified that the synchronized transforma- 
tion is prese nt f or all asymmetries within the previous 
limit in Eq. |lJ). 

Note that a large mixing angle is essential in obtaining 
flavor "equilibration," i.e., that £^ is effectively trans- 
ferred to £f as shown by DHPPRS |lj|. The underlying 
dynamics is simply the adiabatic evolution of the initial 
neutrinos into vacuum mass eigenstates. This would be 
exactly the usual MSW effect 5 if it were not for addi- 
tional complexity of synchronization. We achieve equi- 
libration in the sense that the initial asymmetry is par- 
titioned across the flavors (with the ratio of the final 
and T set by the vacuum mixing angle) . This "equi- 
libration" is simply a MSW transformation that leaves 
the ensemble in a coherent state. This is to be distin- 
guished from equilibration in the conventional sense of a 
completely incoherent or relaxed state, i.e., one produced 
by collisions. 



III. NEUTRINO PROPERTIES AND 
ASYMMETRY TRANSFORMATION 

Since the asymmetry in the electron neutrino number 
is the most stringently constrained, its enhancement due 
to coupling to the other flavors is crucial. The neutrino 
oscillation solution best fitting the observed solar electron 
neutrino flux and spectra is the region of mixing parame- 
ter space named the large mixing angle (LMA) solution, 
with maximum likelihood parameters for two-neutrino 
mixing of order Suiq «4x 10~ 5 eV 2 and sin 2 28 w 0.8. 
Since the LMA mixing is large but not maximal, the first 
mass state (|wi)) is more closely associated with the elec- 



If V T = 0, mb find A sync — <5mj)/132T for £ = 0.05 in agreement 5 Note, however, that for a normal hierarchy we do not have a res- 
with Rcf. |19| . We note that the momentum scale p/T ~ 132 onance — the negative thermal potential makes the u e 's lighter, 

does not determine the character of the solution. and does the same for P e . 
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tron neutrino and the other (\rri2)) less so, and in order 
to enable resonance in the sun, mi < m 2 . 

We also know from atmospheric neutrino observations 
that fi and r neutrino flavors are maximally mixed su- 
perpositions (or nearly so) of two mass states \rri2) and 
\rri3). Therefore, the flavor composition of mass state 
\rri2) is that of a nearly maximal superposition, and com- 
plicates discussion of LMA mixing in neutrino environ- 
ments where the flavor content of \rri2) is of interest. 

However, a powerful simplification can be made given 
maximally mixed (or more generally "similarly coupled" ) 
Ufj, and v T , which allows a linear transformation of the 
3x3 mixing matrix such that one effective flavor state 
\v*) is identically a vacuum mass eigenstate and decou- 
ples from the matter effects [|8) . It is sufficient to follow 
the two remaining states \v e ) and \v*}- This is only jus- 
tified if the momentum state (or more exactly, momen- 



tum distributions) of the superimposed flavors 



and 



\v T ) are indistinguishable. That is, the temperatures and 
chemical potentials of and v r should be equal, i.e., v^, 
v T should be equilibrated. The atmospheric neutrino re- 
sults actually can provide th at v u and v T are equilibrated, 



which we discuss in Section [II B. Strictly speaking, the 
similar-coupling limit is exact only when J7 e 3 = 0, and 
we consider the more general case below. 
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FIG. 4: The evolution of Ji in the synchronized case with 
LMA parameters. As described in the text, the behavior is es- 
sentially a MSW transformation u e <-» v*. The antineutrinos 
I Ji\ evolve identically. The fact that J y is never large demon- 
strates that all of the precession angles are small enough that 
the evolution is dominantly in the x — z plane. The evolution 
of Pi and Pi at the average mom entu m is the same if one 
excludes the self-potentials in Eq. (2.4). 



synchronized mixing angle 



Electron flavor transformation 



26, 



sync 



1 

z 



5m 2 l /2p sync 

V T (j) S ync) 



(3.1) 



Recall that we are interested in whether an asymmetry 
initially present in the poorly constrained or v r will 
convert into a stringently constrained v e jv e asymmetry. 
We can analyze how the LMA mixing parameters evolve 
in the early universe through the effective two-neutrino 
system (z/ e ,i/*). The initial system may be prepared by 
an unspecified leptogenesis mechanism to be in an un- 
mixed state with the asymmetry in v* (£*, 7^ 0) and no 
asymmetry in is e (£* = 0) and remains in this state from 
damping by collisions. We have defined the direction +z 
to correspond to the v e flavor. Taking, for the sake of 
the example, the initial to be negative, the vector I 
will initially point in the +z direction. 

At initially high temperatures, the effective magnetic 
field vector A e ff is dominated by the thermal lepton po- 
tential V T , and is aligned in the — z direction. As the 
universe cools, A e ff rotates away from —z and asymp- 
totes to its vacuum value, which lies close to the +x 
direction (i.e., the angle it makes with the z axis is 20q). 
A large vacuum mixing angle is clearly necessary for this 
MSW transformation to work. 

For an initial asymmetry of |£* » | = 0.05, the evolution 
of the synchronized vector components Ji are shown in 
Fig. |^. The components are driven as a magnetic dipolc 
adiabatically following the evolution of the magnetic field 
A c ff. The evolution of Pi (Pi) at the average mome ntum 
is the same if one excludes the self-potentials in Eq. (2.4). 
The power law growth of J x is simply the evolution of the 



and the growth of J y ~ T results directly from ( 2.11 ) 



Obviously, the transformation occurs when the orienta- 
tion of J rapidly evolves at (A p ) - (V T ) at T ~ 2 MeV. 
The temperature of the transition point in J x scales only 
as (5m 2 ) 1 / 6 , so the results are rather insensitive to the 
uncertainty in the LMA 5m 2 . The antineutrinos evolve 
identically by following the vector — A e fj. The final state 
of the asymmetry after the MSW transformation entering 
the nucleosynthesis epoch is then transferred in propor- 
tion to the vacuum mixing amplitude between the two 
flavors, i.e., J z (l MeV): 



2 I " i! 
1 + cos 26> 



(3.2) 
(3.3) 



(Some care must be taken in interpreting the limits on 
, since the final distributions are not exactly thermal, 
as they are superpositions of Fermi-Dirac distributions 
with different chemical potentials.) Obviously, complete 
"equilibration" or = only occurs for maximal mix- 
ing. The antineutrino chemical potential evolves to the 
values £f = — £| and & — We note that colli- 

sions (which we have neglected) will help make the flavor 
transformation more complete and thus should reduce 
the sensitivity to the mixing angle. 
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B. Mu-tau flavor transformation 

Maximal neutrino mixing as indicated by the Super- 
Kamiokande observations of atmospheric neutrinos has 
nearly identical implications for the evolution of asym- 

Because of the hierarchy 
A p ) ~ (V T ) at the higher tempera- 
ture T ~ 10 MeV. Necessary in driving the flavor evolu- 
tion here is the presence of the remnant thermally pro- 
duced charged muons with energy density 



metries between and v. 

5m itm > KmA. < 



1 



p 2 dp 



P 



1 + cxp 



Jp 2 +ml/T) 



(3.4) 



Though far from the thermal abundance of at T ~ 
20 MeV, real muons remain enough to dictate the flavor 
evolution. Ho wever, since T < m^, the thermal potential 
from Eq. is modified as {E^±) -> j(E^± +p 2 /3E fl ±) 
p9| . We solved the evolution of this case numerically, 
explicitly including the thermal abundance of fi^ 1 , whose 
disappearance accelerates the growth of J x and J v away 
from the power-law growth in the previous LMA case, 
but for simplicity have ignored collisions (Fig. |h. Max- 
imal mixing then gives an equilibration = , which 
allows the application of the simplifying basis transfor- 
mation of the previous section. Inclusion of collisions 
would damp the oscillations at low temperatures but not 
the transformation, as found by DHPPRS |l9[ . 

Interestingly, in the case of evolution without the pres- 
ence of thermal ^ , the evolution is different, with pure 
synchronized vacuum oscillations taking place (rotation 
in the z — y plane) after the Hubble time exceeds the os- 
cillation time. Collisions, which we have omitted, would 
modify the oscillations seen here. 




FIG. 5: The evolution of Ji (Ji are identical) for the mu and 
tau neutrino transformation with and without the inclusion 
of thermal /i pairs. The spiky features indicate real oscilla- 
tions going through zero, and the depth of the spikes on the 
logarithmic scale is an artifact of numerical sampling. Those 
oscillations are real and are determined by the atmospheric 
SrriQ. In the lower panel, J x is zero since the mixing angle is 
maximal. Collisions have been ignored. 



C. Effects of U e3 

The possibility of a nonzero value of U e s obstructs the 
simplifying linear transformation to the basis \v e ) and 
However, nonzero U e s may allow partial equilibra- 
tion of ^,^ T into £ e earlier, at T ~ 5 MeV. For solar 
LMA mixing, significant transformation will always oc- 
cur at T ~ 2 MeV so the value of U e ^ will not alter the 
basic outcome. However, substantial equilibration at 5 
MeV (well before the beta-equilibrium freeze-out) makes 
the general conclusions even more inevitable. 

There are, however, some subtleties associated with 
the sign of Sm^^ — that is, whether the neutrino spec- 
trum has a normal or inverted hierarchy. 6 For a normal 



The sign of the solar 8m 2 is determined by the requirement that 
there be a MSW transition in the Sun, which precludes a reso- 
nance in the early universe (for both neutrinos and antincutri- 
nos). There is, however, no such constraint of the sign of the 
atmospheric 5m 2 . 



hierarchy, the fact the thermal potential makes the i/ e 's 
(and V e 's) lighter implies that that no resonance condi- 
tions can be satisfied in the early universe. With an in- 
verted hierarchy, however, a v e — v* resonance will occur 
when V T ~ Sm 2 tln . We plot in Fig. || level crossing dia- 
grams for neutrinos of the average energy in the absence 
of the self potential. As discussed above, this is a very 
good description of the evolution of the entire neutrino 
distribution. 

The U e 3 mixing angle is constrained to be small [ j30| . 
and as such, coherent evolution will not lead to large fla- 
vor transformation (for the inverted hierarchy, coherent 
evolution through the resonance would swap asymme- 
tries between flavors). However, at T ~ 5 MeV colli- 
sional processes are still highly important and will help 
achieve equilibration. This should be somewhat more 
effective for the inverted case (where the mixing angle 
goes through a maximum) as collisions equilibrate most 
effectively when mixing angles are large. 
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FIG. 6: Level-crossing diagrams for neutrinos of the average 
momentum in the absence of the self-potential. In the upper 
panel we have a normal hierarchy, where the neutrino mass 
eigenstates asymptote to their vacuum values, without ever 
going through a resonance. This is to be contrasted with the 
inverted hierarchy shown in the lower panel where nonzero 
Ue3 leads to a resonance at T ~ 5 MeV. 



IV. NEW CONSTRAINTS 

The electron neutrino asymmetry £ e is limited by its ef- 
fects on the primordial 4 He abundance, Y p . At nucleosyn- 
thesis, nearly all neutrons are incorporated into 4 He nu- 
clei, and Y p production is limited by the neutron fraction, 
set by the freeze-out of beta equilibrium at T ~ 1 MeV. 
The change in the neutron to proton ratio with non-zero 
£ e is simply a Boltzmann factor n/p oc e~£ c ml — £ e . 
And since Y p oc n/p, the uncertainty in the constraint on 
£ e is directly related to the uncertainty in the primordial 
helium abundance AY, 



v 



AYp 

Yn 



(4.1) 



Therefore, one can be very conservative regarding the 
error on the primordial abundance, e.g., AY p « ±0.010 
jn| and still limit A£ e w ±0.04, or equivalently, \L e \ < 
0.03. This method ultimately relies on the uncertainty 
(mostly systematic) in the primordial abundance of 4 Hc. 
Refinement of this constraint may be possible by applying 
CMB priors to BBN predictions combined with reduced 
systematic uncertainties of observed primordial element 
abundances t$2|. 



Analysis of synchronized transformation of neutrino 
asymmetries indirectly translates the constraints on 
to and Q.. In the extreme scenario, an asymmetry 



in (y T ) is equilibrated with v r (y^ 
mixing, such that the state v* has £ 



til. 



) for maximal 
0-5C.r- The 



LMA solution transforms £ M * as Eq. (3.3) so that 



1 — cos 26q 



C ■ 



(4.2) 



For the best-fit LMA mixing angle sin 2 26q w 0.8, the 
limit on an initial asymmetry is ^ T < 0.3. However, 
the LMA mixing angle is not precisely specified. The 
lower end of the 95% confidence level (C.L.) region has 
sin 2 2#o ~ 0.6, for which the limit on the initial asym- 
metry is considerably weaker 7 T < 0.5. The effective 
"2ct" limit therefore is actually an order of magnitude 
larger than that given in DHPPRS since a "small-angle" 
LMA solution reduces the transformation amplitude con- 
siderably. The sensitivity of the KamLAND experiment 
to the LMA parameter space can confirm the LMA pa- 
rameters |$3] and potentially reduce the mixing angle 
uncertainty, and thus improve constraints on the lepton 
number. 



V. DISCUSSION AND CONCLUSIONS 

Due to synchronization by the neutrino self-potential, 
transformation of a large fraction of any asymmetries in 
or v r number to v e is an inescapable consequence of 
the near bimaximal mixing framework emerging for the 
neutrino mass matrix. We have performed a full numer- 
ical integration of the evolution equations in Eq. fl2.4|). 
The numerical solution is nontrivial due to stiff, non- 
linear equations with terms whose time scales vary by 
several orders of magnitude. We confirm the numerical 
results of DHPPRS, and agree that large initial asym- 
metries in and v T are effectively transformed into a v e 
asymmetry, so that the bound from BBN bounds all [|l9| . 

In addition, we have shown numerically that the cou- 
pled evolution of the full-momentum results can also be 
obtained in the average-momentum case when the non- 
linear coupling is neglected. The transfer of neutrino 
asymmetries between flavors occurs identically even when 
i gnor ing the numerically dominant self-potential. In Eq. 
( |2.19| ) we have derived that the self-potential drives a 
synchronization of all momenta to a momentum mode 
p/T = 7T, so that the system by numerical coincidence 
closely follows the average momentum case p/T ~ 3.15. 

We conclude by considering the following implications 
of these results: 



7 Note that we expect this limit would be tighter were we to include 
the effect of collisions. 
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(1) The uncertainty in the lepton number of the uni- 
verse may be reduced by up to two orders of magnitude. 
However, the most conservative limits place 



K 



< 
< 



0.04, 
0.5, 



(5.1) 
(5.2) 



QL e \ < 0.03 and |i„ + L T \ < 0.4). These limits will 
be improved by reducing systematic uncertainties in the 
inferred primordial 4 He abundance and the precise de- 
termination of the baryon density by satellite anisotropy 
experiments Microwave Anisotropy Probe (MAP) and 
Planck |Q. It also may be improved by verification of 
the LMA parameters by KamLAND, particularly if the 
mixing angle is at the large end of the presently allowed 
range. The upcoming data from SNO |35[ will also play 
a very important role in reducing the mixing parameter 
uncertainties. 

(2) Because effectively asymmetries in any neutrino 
flavo r will affect beta equilibrium, the stringent limits 
( |5.l| ) consequentially eliminate the possibility of degen- 
erate BBN |2(J , since an increase the expansion rate with 
large |£ M T | ~ 1 can no longer be compensated by a small 

6 - o.i. 

(3) The above limits on degeneracy in terms of extra 
relativistic degrees of freedom Ai\T„ [see Eq. (1.3)1 are 
impressively small: AN V < 0.004 for the best-fit LMA 
solution, and AN U < 0.2 for the lower limit on the mix- 
ing angle in the LMA solution. DHPPRS suggest that 
AiV„ can be eliminated as a cosmological parameter in 
upcoming fits to the precision CMB data ||. It is cer- 
tainly true that £ can be eliminated, but that is not the 
only possible contribution to AN V . If any nonstandard 
contribution to the relativistic energy density were to be 
detected via the CMB, its origin would be something more 
exotic than degenerate neutrinos, e.g., the decay of a mas- 
sive particle to relativistic species after BBN but before 
CMB decoupling ||. 

(4) It is actually still possible that the upper limit for 



£ e in Eq. (1.4) be fulfilled. Strictly speaking, we have set 
tight new degeneracy limits assuming no non-standard 
contribution to the energy density at the time of BBN. 
It is conceivable that £ e ~ £ M ~ £ T ~ 0.2 if another 
relativistic particle or scalar field contributes the extra 
energy density required to compensate for the large v e 
chemical potential. In this case, flavor-transformation 
improves the current £ MiT limits by at most an order of 
magnitude. Such an unnatural scenario can be detected 
by comparison with the CMB. 

(5) A possible complication to the scenario presented 
here could be mixing with a light sterile neutrino. Ob- 
viously, if the LSND result is confirmed by MiniBooNE, 
then the physics will be much more complicated than as- 
sumed here. If the LSND result is not confirmed, there 
is still the possibility of subdominant mixing to steriles 
that may be difficult to detect in neutrino oscillation ex- 
periments, but which may still play an important role 
in the early universe. Such scenarios have not yet been 
explored. 

(6) A final complication is the yet-unexcluded possi- 
bility of a low reheating temperature (T ~ 1 MeV) p7[ , 
such that the initial conditions of thermal or chemical 
equilibrium for neutrinos for the analysis presented here 
is invalid. Stronger constraints on low-temperature re- 
heating scenarios may be obtained by studying their ef- 
fects on the light element abundances in detail. 
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